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Gravitational waves transport very detailed information on the structure and evo- 
lution of astrophysical sources. For instance a binary system in the early stages of 
its evolution emits a wavetrain at specific frequencies that depend on the charac- 
teristics of the obital motion; as the orbit shrinks and circularize, due to radiation 
reaction effects, the orbiting bodies get closer and tidally interact. This interac- 
tion may result in the excitation of the proper modes of oscillation of the stars, 
and in the emission of gravitational signals that carry information on the mode 
frequencies, and consequently on the equation of state in the stellar interior. These 
phenomena may occur either in solar type stars with orbiting planets and in com- 
pact binaries, and in this lecture we will discuss different approaches that can be 
used to study these processes in the framework of General Relativity. 

In recent years, major progresses have been done in the construction of 
large ground-based interferometers, which will make accessible to observations 
a large frequency region, ranging from a few Hz to a few kHZ: TAMA has 
already performed some observational runs, reaching a sensitivity that would 
enable to detect the coalescence of a neutron star-neutron (NS-NS) binary 
system occuring in our Galaxy with a signal to noise ratio higher than 30, 
and VIRGO, LIGO and GEO600, are almost completely assembled and within 
the end of the year 2002 are expected to start scientific runs El. Resonant 
bar detectors, ALLEGRO, AURIGA, EXPLORER, NAUTILUS, NIOBE, have 
been operational since many years, and are sensitive to a small frequency 
region of about 10-20 Hz, centered at ~1 kHz. As an example, Explorer and 
Nautilus would be able to detect, with signal to noise ratio equal to 1, a burst 
of gravitational waves of amplitude h ~ 4 • 10 -19 , which would correspond 
to a mass of a few units in_10~ 3 M trasforming into gravitational waves 
at the center of our Galaxy 13. And finally, the space-based intereferometer 
LISA, which is expected to fly in about a decade, will enlarge the observational 
window to the low frequency region 10 -4 HzZi/ilO^ 1 Hz. 

Many are the astrophysical sources that, according to the theory of Gen- 
eral Relativity, emit gravitational waves in the frequency region spanned by 
these detectors, and in this review I will discuss the characteristic properties 
of the signals emitted by some of the more interesting sources, and show how 



to compute them. The detectability of a signal depends also on the event rate, 
namely on how many events per year occur in the volume accessible to ob- 
servations, and on the detection rate, which depends on the responce of each 
detector to a specific source; how these rates are determined is, of course, a 
very important issue, but a discussion of these problems is beyond the scope 
of this review, which will be focussed on the information that the gravitational 
signals carry on the structure (either geometrical or internal) of the source, 
on its motion, and on physical processes that may occur in some particular 
situations, like the resonant excitation of the stellar modes. 

In Sec. I, I will shortly describe the quadrupole formalism, which is the 
easiest method to estimate gravitational fluxes and waveforms, and I will dis- 
cuss its domain of applicability. In Sec. II and III, the quadrupole formalism 
will be applied to rotating neutron stars and to binary systems far from co- 
alescence, respectively. In Sec. IV the resonant excitation of the g-modes in 
extrasolar planetary systems will be considered in the framework of a pertur- 
bative approach, and in Sec. V the gravitational emission of coalescing binary 
systems will be studied. 

1 The quadrupole formalism 

The quadrupole formalism is one of the most powerful tools to estimate the 
amount of radiation emitted by a dynamically evolving system. It is based 
on the assumptions that the gravitational field is weak, so that gravitational 
interaction do not dominate, and that the velocities of the bodies involved in 
the problem are much smaller than the speed of light; as a consequence, the 
spacetime metric can be written as the metric of the flat spacetime, n^, plus 
a small perturbation 



and the Einstein equations can be linearized and solved. By a suitable choice 
of the gauge, the perturbation h^ v is found to satisfy a wave equation, and to 
propagate in vacuum at the speed of light; its amplitude is shown to depend 
only on the time variation of the energy density of the source, p(t, x), as follows 
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where q 13 is the quadrupole moment, and the integral extends over the source 
volume. 



An interesting point to stress: the asumption 

v « c — ► flR « c, 

where R is the source extension, implies that the region where the source is 
confined must be small compared to the wavelenght of the emitted radiation, 
Xgw = ^jf- When this constraint is not verified, the quadrupole approach is 
no longer applicable. For instance, the binary system PSR 1931+160, whose 
emission properties will be described in some more detail in the next section, 
is composed of two neutron stars with masses Mi ~ M2 ~ 1.4M Q , moving on 
an eccentric orbit of diameter D ~ 13.8 • 10 10 cm; due to the orbital motion it 
emits gravitational waves at some specific frequencies, and the highest spectral 
line is at vqw ~ 10 -4 Hz. It follows that 

Xgw = ~ 10 14 cm, Xgw >> D, 

vgw 

and the quadrupole constraint is satisfied. But if we want to describe a pul- 
sating neutron star, since the frequency of the fundamental mode of oscillation 
ranges within ~ 2 — 3 kHz, depending on the equation of state (EOS), and 
since the typical diameter of a NS is D ~ 20km, 

c 3 ■ 10 10 cm/ a 7 

^gw = -z — » 1A3 „ — 10 cm, i.e. X G w ~ D, 

vgw 2 — 3 • 10 d Hz 

and the quadrupole formalism cannot be applied. To study these problems, 
alternative and more sophisticated approaches have to be used as we shall see 
in Sec. 0. 



2 Rotating neutron stars 

As a first application of the quadrupole formalism, we shall consider the case 
of rotating neutron stars. Neutron stars can radiate their positive rotational 
energy essentially in two ways. If the star is triaxial, it has a time varying 
quadrupole moment and, according to Eq. (|l|), emits gravitational waves at 
twice the rotation frequency, and with an amplitude which can be parametrized 
in the following wayu 

^4.2xl(r 24 ra 2 (-!— V 1 J ? (<) (2) 
\PJ \10kpcJ 10 45 g cm 2 VICT 6 / w 

where I = 10 45 g cm 2 , is the typical value of the moment of inertia of a 
neutron star, P is the rotation period, and e is the oblateness of the star, here 



given in units of 10~ 6 . These waves could be detected, for instance, by VIRGO 
with one year of integration, if the amplitude of the signal were of the order of 
h r~j 10~ 26 ; this means that neutron stars outside our Galaxy would practically 
be undetectable. In addition, much depends on the value of the oblateness, 
and several studies have tried to set limits on its possible range of variation. 
Gourghoulon and Bonazzola made a first attempt to set constraints on e using 
observational dataB. They considered a number of pulsars, and assuming that 
the observed slowing down of the period is entirely due to the emission of GW, 
gave an estimate of e. They found e [10~ 2 , 10~ 9 ], but of course this can 
only be an upper bound, since we know that rotational energy is dissipated 
in pulsars also by other mechanisms, like the electromagnetic emission and/or 
the acceleration of charged particles in the magnetosphere. Further studies on 
this problem cl established that if the triaxial shape is due to strains in the 
crust of the neutron star, the oblateness could be e ~ 1CP 7 ( TTpj) , where 
a is the strain needed to break the crust-, but its value is quite uncertain; 
for_instance a S [10~ 2 ,10 _1 ] according toQ, and a S [10~ 4 , 10~ 3 ] according 
tou. A time dependent quadrupole moment can also be due to a precession 
of the star's angular velocity around the symmetry axes. In this case the 
radiation is emitted at a frequency v vrec — -^{uj rot + uj prec ) ~ v ro ti but the 
amplitude of the precessional contribution depends on a further parameter, the 
"wobble-angle" between the rotation and the symmetry axes, which is largely 
unknown. From the above discussion, we can conclude that unless astronomical 
observations will set more stringent constraints on the oblateness and on the 
wobble-angle of rotating neutron stars, it is difficult, at present, to say whether 
the radiation emitted by rotating neutron stars will be detectable in the near 
future. 

It is worth mentioning that the detection of gravitational waves depends 
on many cooperating factors; indeed a signal is detectable if its frequency is 
in the range of the detector's bandwith, if its amplitude is high enough to 
be extracted from the noise by a suitable filtering technique, and, last but 
not least, if the computing resources are sufficient to execute this operation. 
For instance, in the case of rotating neutron stars, the estimated number of 
sources m the Galaxy is - 10 9 , of which, ~ 1000 are observed as pulsars, and 
5 are at a distance smaller that 200pc; thus, at least the number of possible 
sources is encouraging; for some of them we have some further information, 
like the location in sky and the spin-down rate. However, we do not know the 
oblateness and since the star is certainly in relative motion with respect to the 
detector, the filters to be used to extract the signal from noise must contain 
the Doppler correction; thus, to detect signals from stars of unknown position 
and characteristics, we should divide the sky in patches as small as possible, 



and search for possible sources by using filters with different parameters, and 
this is extremely expensive! 



3 Binary systems far from coalescence: the quadrupole approach 

The dynamical evolution of a binary system is affected by the emission of 
gravitational waves; due to the energy loss the orbit contracts, the orbital 
velocity increases and the system emits more gravitational energy. The process 
of inspiralling proceeds faster and faster as the system shrinks and finally the 
two bodies coalesce and merge. In this section we shall consider only systems 
that are far from coalescence. Since the two stars are far apart, their internal 
structure is not affected by the tidal interaction, and it is possible to treat 
them as two pointlike masses revolving around their common center of mass. 
In this case the time varying quadrupole moment is entirely due to the orbital 
motion of the two masses, and the emitted waveform is easily calculated by 
the quadrupole formula ([!]). 

In 1975, Hulsc and Taylor applied the quadrupole formalism to. predict 
the slowing down of the period of the binary pulsar PSR 1913+16S. They 
found 4? = —2.4 • 10~ 12 , in excellent agreement with the observed value, 
4^ = —(2.3 ± 0.22) ■ 10 -12 , thus providing the first indirect evidence of the 
existence of gravitational waves. But now the question is: can these waves 
be detected directly? The quadrupole formula shows that when the orbit is 
circular the radiation is emitted at twice the keplerian orbital frequency 
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where M = Mi + Ms is the total mass of the system; if it is eccentric, as in 
the case of PSR 1913+16, waves will be emitted at frequencies multiple of 
V). (in this case, R must be replaced in Eq. (|^) by the semi- major axis), and 
the number of equally spaced spectral lines will increase with the eccentricitytl. 
It is convenient to introduce a characteristic amplitude which can be compared 
to the detectors sensitivity as follows^ 



/'-<»"/.■'•> = \/§ (h$) 2 + (h%!) 2 • (4! 
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a PSR 1913+16 is composed of two very compact stars, with masses Mi = 1.4411 Mq 
and M2 = 1.3874 Mq, revolving around their center of mass with an eccentric orbit 
(e = 0.617139), and keplerian frequency = (u}] s /2w) = 3.583 ■ 10~ 5 Hz. The system is at 
a distance D = 5 kpc from Earth. 
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Figure 1: The characteristic amplitude of the wave emitted by PSR 1913+16 is plotted 
versus the harmonic index n. The spectral lines are emitted at frequencies multiple of v^. 

where (fig^) 2 and (h 1 ^) 2 , are the square of the n— th Fourier component of 
the two independent polarizations, averaged over the solid angle 

& 2 = ^J dn \hee{nv k ,r,6,cl>)\\ (h<g) 2 = ±- J 'dfi \ h H {nv k , r, 6, 0)| 2 , 

and the factor •v/2/3 takes into account the average over orientation. In the 

previous expressions, the waveforms have been expanded in Fourier series 

Ht) = EZ-Z h(nu; k )e- 2in ^/ p , and h(nu k ) = £ J Q P h(t)e+ 2m ^/ p dt, where 

LO k = 2lTV k . 

The characteristic amplitude h c emitted by PSR 1913+16 is shown 
in figure 1 as a function of the harmonic index n. The maximum occurs 
at v m ax — 1-44 • 10 -4 Hz, and the amplitude of the corresponding line is 
h - If) -23 

llj c max ±yJ 

Thus the emission of this system is in the bandwidth of the space-based in- 
terferometer LISA, but the expected sensitivity curve of this instrument shows 
that the signal is too low to be detectable even with five years of integration. 

However, other interesting sources of radiation of this kind exist in our 
Galaxy. For instance, in a recent paper the gravitational emission of cata- 
clysmic variables has been considered Ej. These systems are semi-detached 
binaries of low mass and very short period, in which the primary star is an 
accreting degenerate white dwarf, and the secondary is usually a late type-star 
filling its Roche lobe and transfering matter on the companion. In table I the 



Table 1: Cataclysmic variables with component masses mi and rri2, gravitational wave 
frequency vgw, and gravitational wave amplitude h c . 



name 


mi(M ) 


m 2 (M ) 


vgw 


(Hz) 


log /l c 


U Gem 


1.26 


0.57 


1.31 • 


10" 


4 


-20.8 


IP Peg 


1.15 


0.67 


1.46- 


io- 


-4 


-20.9 


HU Aqr 


0.95 


0.15 


2.67- 


io- 


-4 


-21.3 


VW Hyi 


0.63 


0.11 


3.12- 


io- 


-4 


-21.3 


EX Hya 


0.78 


0.13 


3.39- 


io- 


-4 


-21.4 


WZ Sge 


0.45 


0.058 


4.08- 


io- 


-4 


-22.1 


ST LMi 


0.76 


0.17 


2.93- 


io- 


-4 


-21.4 


SW UMa 


0.71 


0.10 


4.93- 


io- 


-4 


-21.6 


Z Cha 


0.84 


0.125 


3.11 • 


io- 


-4 


-21.5 


V 436 Cen 


0.7 


0.17 


3.70- 


10" 


-4 


-21.6 



masses, emission frequencies and characteristic wave amplitude are listed for 
some of these short period systems; among the systems analized in ref. 11, 
we have selected only those that emit gravitational waves large enough to be 
detectable by LISA with one year of integration. 

4 Extrasolar planetary systems 

In recent years a large number of extrasolar planetary systems have been dis- 
covered at very short distance from usli3, which exhibit some very interesting 
properties. They are very close to Earth (less than 10-20 pc), and they are 
formed by a solar type star and one or more orbiting companions, which can 
be planets, super-planets or, in some cases, brown dwarfs. Many of these com- 
panions are orbiting the main star at such short distance that conflicts with 
the predictions of the standard theories of planetary formation and evolution. 
For instance, some planets [] have orbital periods shorter than three days (for 
comparison, Mercury's period is 88 days). In the light of these findings, it is 
interesting to ask the following question: is it possible that a planet orbit the 
main star so close as to excite one of its modes of oscillations? How much 
energy would be emitted in the form of gravitational waves by a system in 
this resonant condition, compared to the energy emitted because of the orbital 
motion? For how long would this condition persist? These questions have been 
analysed in two recent paper Lt3, whose results I shall briefly summarize. First 



of all we computed the frequencies of pulsation of the modes of a solar type 



''Here and in the following we shall indicate as "planet" also brown dwarfs and super planets. 



Table 2: Values of the frequencies (in fiHz) of the fundamental and of the first g-modes of 
oscillation of a solar-mass polytropic star with n = 3, for 1 = 2. 



f-mode ffi ff 2 ff 3 ff 4 ff 5 gg g 7 fls ff9 ffio 
Vi 285 221 168 135 113 97 85 75 68 62 57 



star, by integrating the equation of stellar perturbations in the framework of 
General Relativity. We considered a polytropic model of star, p = Ke 1+1 ^ n 
with n = 3, adiabatic exponent 7 = 5/3, central density eo = 76 g/cm 3 , and 
c 2 eo/po p,5.53 • 10 5 . This model gives a star with the same mass and radius as 
the Sun El In table 2 the values of the frequencies of the fundamental mode 
and of the first g-modes of the star are tabulated for £ = 2. We do not include 
the frequencies of the p-modes because they are irrelevant to the following 
discussion. 

If a planet moves on an orbit of radius R (we shall assume for sim- 
plicity that the orbit is circular), the keplerian orbital frequency is — 

_L 2l M^+ M p) ; w here M* and M p are the mass of the star and of the planet, 
respectively. The planet can excite a mode of frequency v>i, only if vi is a 
specific multiple of the orbital frequency, i.e. 

Vi=£vk, (5) 
where I is the considered multipolc. Introducing the dimensionless frequency 
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» R * ' 

where i?* is the radius of the star, the resonant condition (||) can be written 



as 



n 2 / a/, m 1/3 



1 



R*. (6) 



Thus, Ri is the value of the orbital radius for a given mode to be excited; for 
instance, using Eq. (||) we find that in order to excite the fundamental mode 
of the considered star, the planet should move on an orbit of radius smaller 
than 2 stellar radii, therefore we first need to verify whether a planet can move 
on such close orbit without being disrupted by the tidal interaction. This is 
equivalent to establish at which distance the star starts to accrete matter from 



the planet (or viceversa), i.e. when the planet or the star overflow their Roche 
lobes |. 

It is an easy exercise to compute the radius of the Roche lobe, Rrl, in 
newtonian theory, which fully applies to the case we are examining. If we 
choose a reference frame co-rotating with the two masses, with origin on the 
center of mass and such that the x-y plane coincides with the orbital plane, 
and the x-axis joins the two masses, the newtonian potential of the star-planet 
system can be written as 

TTt \ GM P GM * 1 2, ,2 

U(x,y) = p - - -^|x| 2 , 



where = (^ M M + M R,0j , x p = ^— M A ^* M R, 0^ , R is the orbital radius of 
the planet, and the last term is the centrifugal potential. A gas particle at a 
large radius on the planet may feel a gravitational force from the star that is 
comparable to, or even larger, than that from the planet. In that case, the gas 
particle is unstable against being transferred from the planet to the star. This 
can be understood by plotting contours of equal gravitational potential: if the 
potential is large and negative, the equipotential surfaces are two spheroids, 
one centered on M p , the other on M*. As the potential increases the spheroids 
deform, and finally touch in one point lying on the axis which joins the star 
and the planet; the surface they form in this configuration is the Roche lobe 
surface, the unique point where the potential energy contour intersects itself is 
called the inner Lagrange point, and the curve which encloses the Roche lobes 
is said "first Lagrangian curve" . If the planet expands to fill its Roche lobe, the 
matter near the inner Lagrange point will spill over onto the star. Similarly, if 
the star fills its own Roche lobe it can accrete matter onto the planet. We can 
define the Roche lobe radius for the planet, Rrl, as the radius of the circle 
centered on the planet and tangent to the first Lagrangian curve (similarly for 
the star). Rrl can be found numerically, and it depends exclusively on the 
planet orbital radius R and on the star and planet masses. 

The procedure to establish if a planet is allowed to orbit at a given distance 
from the star exciting its modes is the following: 

- we assign a value of R — Ri which corresponds to the excitation of a given 
mode, according to Eq. (||) ; 

- we determine the Roche lobe radius Rrl{Ri 1 M p , M*) and compute the di- 

c Following the analysis in ref. Q, temperature effects, which may provoque the melting of 
the planet or the evaporation of its atmosphere, can be shown to be less stringent than the 
Roche-lobe limit. 



Table 3: The critical ratio PRL is given for three planets with mass equal to that of the 
Earth (Me), of Jupiter (Mj) and of a brown dwarf with Mbd = 40 Mj, and for the different 
modes. This critical ratio corresponds to the minimum mean density that a planet should 
have in order to be allowed to move on an orbit which corresponds to the excitation of a 
g-mode, without being disrupted by tidal interactions. 



Prl 





9i 


92 


93 


34 


35 


36 


97 


98 


99 


3io 


M E 


12.5 


7.21 


4.64 


3.24 


2.38 


1.83 


1.45 


1.17 


0.97 


0.82 


Mj 


13.0 


7.49 


4.82 


3.37 


2.48 


1.90 


1.51 


1.22 


1.01 


0.85 


Mbd 








3.68 


2.71 


2.08 


1.65 


1.34 


1.11 


0.93 



mensionless quantity 



Rrl given by 

Rrl = Rrl/ Ri', 



(7) 



- if the planet radius would be equal to Rrl, its average density would have a 
critical value, prl, given by 



Prl 



h&RL 



(8) 



Thus, if p p > prl the planet will not fill its Roche lobe and will not accrete 
matter on the star. Using Eqs. Eqs. (H),(@) and (Q), the value of the critical 
density can be rewritten as 



M r , 



Prl 



InRl 



r^2 
ft? 



Mp 



(9) 



R 



RL 



or, normalizing to the mass of the central star, 



Prl 
P* 



M p /M* 



M p /M+) R RL ' 



P{1 



(10) 



In conclusion, a planet can excite the ith-mode of the star without overflowing 
its Roche lobe, only if the ratio between its mean density and that of the 
central star is larger than the critical ratio (|10| ). 

We have computed the critical ratio (|10|) assuming that the central star 
has three different companion: two planets with the mass of the Earth and of 
Jupiter, and a brown dwarf of 40 jovian masses, Me, Mj and Mbd, respec- 
tively, imposing that they are on an orbit which corresponds to the resonant 



excitation of a given mode of the solar type star. We do not consider smaller 
planets because they produce gravitational signals that are too small to be in- 
teresting. The results of this analysis are shown in Table [| where we tabulate 
the value of for the three considered companions. 

A planet like the Earth has a mean density such that PEarth/ Po = 3.9, 
whereas for a planet like Jupiter Pj/pq — 0.9; thus, from Table || we see that 
an Earth-like planet can orbit sufficiently close to the star to excite g-modes 
of order higher or equal to n = 4, whereas a Jupiter-like planet can excite 
only the mode <?io or higher. We find that in no case the fundamental mode 
can be excited without the planet being disrupted by accretion._Accordmg to 
the brown dwarf model ("model G") by Burrows and Liebertlij an evolved, 
40 Mj brown dwarf has a radius Rbd = 5.9 • 10 4 km, and a corresponding 
mean density pbd = 88 g-cm~ 3 ; consequently pbd/pq = 64, and this value is 
high enough to allow a brown dwarf companion to excite all the g-modes of 
the central star. However, we also need to take into account the destabilizing 
mechanism of mass accretion from the central star. This imposes a further 
constraint, and this is the reason why the slots corresponding to the excitation 
of the g-modes lower that g4 in the last row of Table | are empty. 

Having established that some g-modes of the central star can, in principle, 
be excited, we turn to the next question, i.e. : how much energy would be 
emitted in gravitational waves by a system in this resonant condition, com- 
pared to the energy emitted because of the orbital motion? The amplitude 
of the wave emitted because of the orbital motion can be computed by us- 
ing the quadrupole formula (Q); in the case of circular orbit the characteristic 
amplitude (^) becomes 

h Q (R l ) = ^ M^R,) 2 . (11) 
r V 15 

The values of this amplitude for the three planets considered above, are shown 
in table ||. However, if the stellar modes are excited the star will pulsate 
and emit more gravitational energy than that which can be computed by the 
quadrupole formula. In order to evaluate the resonant contribution we need to 
follow a different approach which takes into account the internal structure of 
the main star. One possibility is to use a perturbative approach which consists 
in the following. The main star is assumed to be an extended body, whose equi- 
librium structure is described by an exact solution of the relativistic equations 
of hydrostatic equilibrium; the planet is considered as a pointlike mass which 
induces a perturbation on the gravitational field and on the thermodynamical 
structure of the star. Under these assumptions we solve the equations of stel- 
lar perturbations to compute the characteristic amplitude of the gravitational 



Table 4: The amplitude of the gravitational signal emitted when the three companions 
considered in Table H move on a circular orbit of radius Ri, such that the condition of 
resonant excitation of a g-mode is satisfied, is computed by the quadrupole formalism (see 
Eq. for the modes allowed by the Roche-lobe analysis. The planetary systems are 

assumed to be at a distance of 10 pc from Earth. 
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ve emitted when the planet moves on a circular orbit close to a resonance 
In the following, we shall call this amplitude hf>, to indicate that it has 
been computed by solving the relativistic equations of stellar perturbations. 

We have computed Hr for I = 2 (which is the relevant contribution), 
assuming that the planet moves on a circular orbit of radius Rq — Ri + AR, 
where Ri is the radius corresponding to the resonant excitation of the mode Qi, 
for the modes allowed by the Roche lobe analysis. We find that, as the planet 
approaches the resonant orbit, Iir grows very sharply. In figure ^| we plot the 
ratio hfj(R ) /hq{Ri) as a function of AR, to see how much the amplitude of 
the emitted wave increases, due to the excitation of a g-mode, with respect to 
the quadrupole emission. The plot is done for the modes g4, g7, and g 10 , and 
it shows a power-law behaviour nearly independent of the order of the mode. 
It should be stressed that hn(Ro)/hQ(Ri) is independent of the mass of the 
planet, but depends on the selected stellar model. From figure || we see that as 
the planet approaches a resonant orbit the amplitude of the emitted wave may 
become significantly higher than that computed by the quadrupole formalism, 
Hq. Thus the next question to answer is how long can a planet orbit near 
a resonance, before radiation reaction effects move it off. Indeed, the loss of 
energy in gravitational waves causes a shrinking of the orbit of a planetary 
system, and the efficiency of this process increases as the planet approaches a 
resonant orbit. We shall now compute the time a planet takes to move from an 
orbit of radius Rq = Ri + AR, where the amplification factor hn/hq has some 
assigned value, to the resonant orbit Ri, because of radiation reaction effects. 
This timescale will indicate whether a planet can stay in the resonant region 
long enough to be possibly observed. On the assumption that the timescale 
over which the orbital radius evolves is much longer than the orbital period 
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Figure 2: The logarithm of the ratio h l R 2 (Ri + AR)/hQ(Ri) is plotted as a function of A_R 
for the modes g4, g7, and gio. 



(adiabatic approximation), the orbital shrinking can be computed from the 
energy conservation law 



(12) 



where E is the energy per unit mass of the planet which moves on the geodesic 
of radius Rq 

-1/2 

(13) 



Rq 



Ro 



and (^^p 1 ) is the energy emitted in gravitational weqzes, computed numer- 
ically by using the perturbative approach (see ref. 13 for details). Since 
<f > = <^> / <^>> using Eq. © Eq. © gives 



dR Q 
~~dT 



2R 2 (1 - 3M+/R Q f/ 2 I dE GW \ 
M P M* (l-6M*/iJo) \ dt I 



(14) 



from which the time needed for the planet to reach the resonant orbit can be 
computed 



AT = 



M P M* 



Ri (1 - QA-hlR ) dR 
fl i+ AH(l-3M 4 / J R ) 3 / 2 ^ (^fP)' 



(15) 



It should be noted that since ( dE ^ t w ) is proportional to M£, AT is longer 
for smaller planets. We have computed AT for the three companions Me, Mj 
and Mbd, f] by the following steps: 

- we assign a value of the orbital radius R — R t which corresponds to the 
excitation of a g-mode allowed by the Roche lobe analysis 

- we compute the radius of the orbit Rq = Ri + AR for which the amplification 
factor A = hn(Ro)/hQ(Ri) has an assigned value 

- we compute the energy radiated in gravitational waves, ( dE ^ t w ) i as a function 
of the orbital radius, in the range between Rq and Ri 

- using Eq. (|l^) we compute AT, which tells us how long can the planet 
orbit the star in the resonant region between Rq and Ri, emitting a wave of 
amplitude higher than [A x hq^Ri)]. 

The results are summarized in table ||. These data have to be used together 
with those in table || as follows. Consider for instance a planet like the Earth, 
orbiting its sun on an orbit resonant with the mode g4. According to the 
quadrupole formalism it would emit a signal of amplitude hQ^R^) = 3 • 10 -26 , 
(table [|, first row) at a frequency vow = 1.13 • 10~ 4 Hz (table || first row) 
. The data of table |5|, which include the resonant contribution to the emitted 
radiation, indicate that before reaching that resonant orbit of radius R4, the 
Earth-like planet would orbit in a region of thickness AR — 1.7 km slowly 
spiralling in, emitting waves with amplitude hji > 10/iQ(i?4) = 3 • 10~ 25 , for a 
time interval of 5.4- 10 6 years, and that while spanning the smaller radial region 
AR = 312 m, the emitted wave would reach an amplitude h,R > 50/iq(-Z?4) = 
1.5 • 10~ 24 , for a time interval of 3.8 • 10 4 years. A jovian planet, on the 
other hand, could only excite modes of order n — 10 or higher, which would 
correspond to a resonant frequency of vgw — 5.7 • 10 -5 Hz and a gravitational 
wave amplitude greater than 6 • 10 -23 for 310 years (AR = 8 m), or greater 
than 3 • 10 -22 for ~ 2 years (AR = lm). From these data we see that the 
higher the order of the mode, the more difficult it is to excite it, because the 
region where the resonant effects become significant gets narrower and the 
planet transits through it for a shorter time. Much more interesting are the 
data for a brown dwarf companion. In this case, for instance, the region that 
would correspond to the resonant excitation of the mode g4 (AR = 1.7 km), 
with a wave amplitude greater than 3.9 • 10 -21 , would be spanned in 430 years, 
whereas the emitted wave would have an amplitude greater than ~ 2 • 10 
(AR = 312 m) over a time interval of ~ 3 years. The results of this study 
show that if a brown dwarf would orbit a solar type star at such distance that 



d The data in table 5 have been obtained by repeating the calculations described in pef. lij, 
after correcting an error found in the numerical code. The qualitative results of ref. |]j arc, 
however, correct. 



Table 5: For each mode allowed by the Roche lobe analysis (see text), we give the frequency 
of the wave emitted when a planet moves on an orbit resonant with a g-mode of its sun 
(column 2). When the planet spans a radial region of thickness AR (column 3) approaching 
a resonant orbit, due to the stellar pulsations the amplitude of the emitted wave is amplified, 
with respect to the quadrupole amplitude, by a factor greater than A (column 4). In the 
last three columns we give the time interval AT needed for the three companions to span 
the region AR and reach the resonance. 



Mode 


v G w (A*Hz) 


AR{m) 


A 


A?Myrs) 


AT BD (yrs) 


ATj(yrs) 


g4 


113 


1700 


10 


5.4 


10 B 


4.3- 


10 2 








312 


50 


3.8 


10 4 


3.0 






g5 


97 


616 


10 


2.7 


10 B 


2.1 • 


10 2 








113 


50 


1.9 


10 4 


1.5 






g6 


85 


240 


10 


1.4 


10 6 


1.1 • 


10 2 








44 


50 


9.6 


10 3 


7.5 • 


io- 1 




g7 


75 


98 


10 


7.0 


10 5 


55 










18 


50 


5.0 


10 3 


3.9- 


io- 1 




gs 


68 


42 


10 


3.8 


10 5 


30 










8 


50 


2.5 


10 3 


1.9- 


io- 1 




g9 


62 


18 


10 


1.9 


10 5 


15 










3 


50 


1.5 


10 3 


1.2 • 


io- 1 




gio 


57 


8 


10 


1.0 


10 5 


7.8 




3.1 • 10 2 






1 


50 


6.7 


10 2 


5.3- 


IO" 2 


2.1 



a g-mode is excited, the emitted radiation may be strong enough, and for a 
sufficiently long time interval, to be detectable by LISA. 

5 Coalescing binary systems 

The main target of ground-based interferometers is to detect the signal emitted 
during the late phases of inspiraling of a binary system. Assuming that the two 
coalescing bodies are two point masses, Mi and M2, by using the quadrupole 
formalism and including the radiation reaction effects (along the lines discussed 
in Sec. 4), it is possible to show that the loss of gravitational energy induces 
a circularization of the orbit t3, and that the radius decreases according to 

1/4 ^ -R 4 

the law R(t) = R in (1 - t/t coa i) ; where t coal = ^jffi is the time of 
the final coalescence, and [i — M\MijM is the reduced mass of the system 
of total mass M — Mi + M2. The orbital frequency consequently increases, 
and so does the frequency of the emitted wave, with a time dependence given 



. r , . , ]3/8 

by v — — 2§6 ^M 2 /a Jt — T 1 *) • Thus, the gravitational signal emitted by 

a coalescing system resembles the chirp of a singing bird. 

At some point of the evolution, for instance ~ 10 s years after formation 

if the binary system is composed of two neutron stars, the frequency of the 

emitted signal enters the bandwidth of the ground based interferometers, and in 

about 15-20 minutes (depending on the mass of the stars) sweeps the frequency 

region ranging from ~ 10Hz up to ~ 900Hz. The amplitudes of the two 

polarizations of the wave emitted during this fast inspiralling phase are 

h+ = ^ '- cos(2tt^), h x = ±— i 1 S in(27rz/t), 

r r 

(16) 

where i is the angle of inclination of the orbit to the line of sight, and the 
emitted energy per unit frequency is given by = 2L j-^ A ^i/3 ■ 

This description of the coalescence based on the quadrupole formalism 
breaks down when the two stars are so close that tidal effects become impor- 
tant. For instance, for binaries with one white dwarf companion, the max- 
imum frequency, i.e. the minimum distance between the two stellar compo- 
nents, is set in order to cut-off the Roche-lobe contact stage. In fact, the 
mass transfer from one component to its companion transforms the original 
detached binary into a semi-detached binary. This process can be accompa- 
nied by loss of angular momentum with mass loss from the system and the 
above description cannot be applied. In this case, the mimimum orbital sepa- 
ration that two white dwarfs can reach is given by r w d(Mi) +r w d(^2)i where 

r wd (M) = 0.012i? Q y / (M/1.44M Q )" 2/3 - (M/1.44M ) 2/3 is the approximate 
formula for the radius of a white dwarf0. If the system is composed, say, by 
two white dwarfs of 0.9 Mq, the cutoff frequency is vgw max ~ 0.1 Hz. These 
signals are therefore candidates for detection by spaceborne interferometers. 

The situation is different for NS-binaries or BH-binaries because, being the 
two bodies extremely compact, the pointlike approximation described above 
works remarkably well up to much higher frequencies, where the post-newtonian 
formalism can be applied to correct the waveforms emitted during the latest 
phases of the inspiral; however, even this more refined waveforms fail to cor- 
rectly reproduce the emitted signal when the innermost stable circular orbit 
(ISCO) is approached, particularly when the coalescing bodies are neutron 
stars. Indeed, in a recent paper we have shown that due to a resonant exci- 
tation of the modes of oscillation of the stars, the emitted energy may be en- 
hanced with respect to the orbital contributionlla. This phenomenon emerges 
only if we use an approach in which the stars are allowed to have an internal 



Table 6: Parameters of the polytropic stars we consider in our analysis: the polytropic 
index n, the central density, the ratio ao = c 2 eo/po of central density to central pressure, 
the mass, the radius and the ratio R/M (in geometric units). The central density is chosen 
in such a way that the stellar mass is equal to IAMq, except for model A, the mass of which 
is about one solar mass. 



Model number 


n 


e (g/cm 3 ) 


a 


M (Mq) 


R (km) 


R/M 


A 


1.5 


1.00 x 10 15 


13.552 


0.945 


14.07 


10.08 


B 


1 


6.584 x 10 14 


9.669 


1.4 


15.00 


7.26 


C 


1.5 


1.260 x 10 15 


8.205 


1.4 


15.00 


7.26 


D 


1 


2.455 x 10 15 


4.490 


1.4 


9.80 


4.74 


E 


1.5 


8.156 x 10 15 


2.146 


1.4 


9.00 


4.35 



structure, which plays a dynamical role by exchanging energy with the gravi- 
tational field. For instance, inlia we have used the same perturbative approach 
described in Sec. || to study the emission of extrasolar planetary systems: we 
assume that one of the two neutron is a "true" neutron star, whose equilib- 
rium structure is governed by the TOV equations of hydrostatic equilibrium, 
whereas the second star is a point mass which perturbes the gravitational field 
and the fluid of the extended companion. To extract the effects that the inter- 
nal structure of the star has on the gravitational emission we have selected a set 
of parameters which allow to encompass a reasonable range of stellar models, 
and to cover most of the range of structural properties obtained with realistic 
EOS. The values of the parameters are given in Table ^. For each model, we 
have integrated the equations of stellar perturbations in General Relativity, 
and we have computed the energy flux emitted in gravitational waves, E R , 
assuming that the perturbing point mass, mo, is moving on a circular orbit of 
radius Rq, with orbital velocity v, and semilatus rectum p, given by 

v = (M Wfc ) 1/3 = 4=, (17) 

where M is the mass of the star. The results are shown in Fig. |^, where 
we plot the normalized energy flux, P = E R /E N , where E is the orbital 
contribution computed by the quadrupole formalism 

E« - § . (18, 

P(v) has been obtained by adding the contributions of different I and m, with 
2 < I < 7, and it is plotted in Fig. || as a function of the orbital velocity, 
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Figure 3: The normalized energy flux, P(v), is plotted as a function of the orbital velocity 
for the stellar models given in Table and for a black hole. For model D and for the black 
hole the curves extend up to the velocity which correspond to the ISCO, whereas for the 
other models they stop when the mass mo reaches the surface of the star. The sharp peaks 
indicate that, for different values of the harmonic index i, the fundamental quasi-normal 
modes of the star are excited if the orbital frequency satisfies the resonant condition (H) ; the 
curve of the stellar model B has a peak at high v which correspond to the excitation of the 
first p-mode for 1 = 2. The most compact model E is not shown in the figure because at this 
scale it is indistinguishable from the black hole. 



for the models of star we have considered; for comparison, we also plot P{v) 
computed in the case when the perturbed object is a black hole. In this case 
the plot extends up to the velocity which corresponds to the ISCO, Rq = 6M, 
whereas for the stellar models it stops when mo reaches the surface of the star. 
From Fig. ||we see that sharp peaks appear if the central object is a star: they 
correspond to the excitation of the fundamental quasi- normal modes of the star 
for different values of the harmonic index £. In the case of model B the first 
p-mode for I = 2 is also excited. Since the peaks of the mode excitation are 
very high, the scale chosen on the vertical axis of Fig. || makes the response of 
the black hole to appear as a flat line. The reason is that, since the frequency 
of the lowest quasi-normal modes of a black hole are much higher than those 
of a star with the same mass, the circular orbit that would excite them would 
have a radius smaller than 6Mbh- Thus, in the range of v considered in Fig. 
| the energy flux emitted by the black hole is due essentially to the orbital 
motion. In Fig. ||, we show a zoom of Fig. || restricted to the region v < 0.28, 
which is far enough from the resonant orbits (except that for model A). In this 
case we can appreciate the differences between the emission of different stellar 
models and that of a black hole. If the orbital velocity is smaller than 0.16 
all curves are practically indistinguishable. Fig. || shows that the normalized 
energy fluxes emitted by different stellar models have a different slope, and 
are always larger than the flux emitted by the black hole. The increase in the 
energy output at orbital velocities larger than v = 0.16 is just an effect of the 
resonance tail. 

6 Concluding Remarks 

In this review we have shown how to estimate the characteristic features of the 
gravitational signals emitted by some interesting astrophysical sources. 

Chances to detect these waves crucially depend on the detailed knowl- 
edge of the waveform; indeed the standard technique used to extract a non- 
continuous signal from the noisy data of a detector is the matched filtering 
technique, whose performances are very sensitive to a mismatching of the pa- 
rameters: for this reason the templates must be as accurate as possible. Since 
the most interesting signals are emitted in spacetime regions where the non 
linearity of the gravitational interaction plays an important role, the ultimate 
template would be that obtained by integrating the fully non linear equations 
of General Relativity. However, fully nonlinear simulations of phenomena like 
the gravitational collapse or the coalescence of binary systems are, at present, 
at a preliminary stage, although major progresses in the field are underway. 
Thus, for the time being we need to rely on the results obtained by approxi- 




Figure 4: The normalized energy flux, P(v), is plotted as in Fig. gj, but for a smaller orbital 
velocity range, such that the peaks due to the excitation of the stellar modes are excluded. 



mate methods, like the quadrupole formalisms and the perturbative approach, 
to mention those we have described in this paper. There are many things that 
we learn from these approaches. Consider for instance the evolution of a binary 
system: the quadrupole formalism, applied to point masses, allows to deter- 
mine with high accuracy the gravitational signal emitted when the two bodies 
are far apart, the frequencies and the expected amplitudes. Including radiation 
reaction effects and using the adiabatic approximation, we can predict how the 
orbit evolves, and how much time it takes to become circular. However, if the 
system is composed by a solar type star and a close planet (or a brown dwarf 
companion), new phenomena may occur, like the resonant excitation of the 
g-modes of the star, which can be accounted for only if the internal structure 
of the star is included in the picture; in this case, the perturbative approach 
provides the tools needed to evaluate how much energy is radiated in gravita- 
tional waves because of the stellar oscillations, at which frequencies, and for 
howlong can a system be in a resonant condition. 

If the binary system is formed by very compact objects, like neutron stars 
or black holes, the quadrupole approach works up to much higher frequencies, 
since in this case the masses can really be treated, to a large extent, as point 
masses. During the latest phases of the inspiralling, let say during the 10- 
20 cicles that preceed the merging when the signal frequency is of the order 
of 900-1000 Hz, the non linearity of the interaction produces deviations from 
the quadrupole signal which can be estimated_by introducing post-newtonian 
corrections at higher order of (v/c) (see e.g. Ej and references quoted there). 
In the standard post-newtonian (PN) approach, the coalescing bodies are still 
considered as point masses, and this is good enough for black holes. However, 
in the case of NS-binaries the fundamental mode of the stars may be excited in 
the same region, and this excitation may produce corrections to the quadrupole 
signal that are comparable or higher than the PN corrections, and depend on 
the internal structure of the stars; as we have shown in Sec. [|, these phenomena 
can be studied again by a perturbative approach, which has been shown to work 
remarkably Jisaell, even when the two stars are at a distance as short as three 
stellar radii E£l. 

An accurate detection of the signal emitted during the last few cicles before 
coalescence may enlighten unresolved issues related to the internal composition 
of neutron stars and to the state of matter at supernuclear density. These 
effects, however will not be seen by the first generation of interferometers, 
for which the most likely source to be detected is the coalescence of BH-BH 
binaries with masses of the order of 20 — 30 M EHl. But there is hope for 
the future: advanced detectors are beeing proposed like EURO, a european 
project now under a feasibility study, that should be extremely sensitive at 



high frequency, and would allow to study the processes we have described. 
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